Typeset with jpsj2.cls <ver.l.2> 



Letter 



in 
o 
o 

o 

Q 



X3 
O 

> 

ON 

cn 



in 
o 



O 

o 



X 



Coherent Behavior and Nonmagnetic Impurity Effects 
of the Spin Disordered State in NiGa2S4 
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Nonmagnetic impurity effects of tfie spin disordered state in the triangular antiferromag- 
net NiGa2S4 was studied tlirough magnetic and thermal measurements for Nii_a;ZniGa2S4 
(0.0 < a; < 0.3). Only 1 % substitution is enough to strongly suppress the coherence observed 
in the spin disordered state. However, the suppression is not complete and the robust feature 
of the dependent specific heat and its scaling behavior with the Weiss temperature indicate 
the existence of a coherent Nambu-Goldstone mode. Absence of either conventional magnetic 
order or bulk spin freezing suggests a novel symmetry breaking of the ground state. 

KEYWORDS: triangular lattice, NiGa2S4, spin disordered state, Nambu-Goldstone mode 



Geometrically frustrated magnets have attracted in- 
terest because of the possible emergence of novel mag- 
netic phases at low temperatures by suppressing con- 
ventional magnetic order. Among such magnets, one of 
the simplest and most fundamental structures is the 
two-dimensional (2D) triangular lattice that has a sin- 
gle magnetic ion per unit cell. The 2D triangular lat- 
tice antiferromagnets (AFMs) have long attracted in- 
terest since a quantum spin disordered state was the- 
oretically proposed more than three decades ago.^ While 
it is believed that the triangular AFM with nearest- 
neighbor coupling exhibits the so-called 120° order, ^""^ 
recent theories suggest that the exchange interaction be- 
yond nearest-neighbor such as longer range and multiple- 
spin exchange interactions may lead to a quantum dis- 
ordered state. For experiment, only a few quasi-2D 
triangular AFMs with low spin, i.e. S < 1, have been 
reported, such as a solid '^He thin film'' and an organic 
material with a distorted triangular lattice.^ 

Recently, a new quasi-2D S = 1 triangular AFM 
NiGa2S4 has been discovered as the first example of a 
bulk low-spin AFM with an exact triangular lattice.^ In- 
terestingly, despite antiferromagnetic (AF) coupling of 
~ 80 K, neither long-range order nor bulk spin freezing 
has been detected down to 0.35 K. Instead, magnetic, 
thermal and neutron diffraction measurements indicate 
the formation of a gapless spin disordered state below 
about 10 K. The quadratic temperature dependence of 
the specific heat, as well as the temperature independent 
behavior of the susceptibility at T — * reveal coherence 
of a gapless linearly dispersive mode in two dimensions. 
Such a spin disordered state is presumably sensitive to 
perturbation such as impurity doping. In order to clar- 
ify the mechanism of the formation of the unusual low 
temperature state, it is highly important to examine ex- 
perimentally the stability against impurities of the spin 
disordered state and its gapless linearly dispersive mode. 

In this letter, we report the nonmagnetic impurity ef- 
fects of NiGa2S4 based on the magnetic and thermal 
measurements for the Zn substituted insulating materi- 
als Nii_a;Zna;Ga2S4. Only 1 % is enough to substantially 



suppress the coherence of the spin disordered state in 
NiGa2S4. However, the Zn substitution never completely 
suppresses the coherence, but induces a "defect" compo- 
nent of weakly coupled S — 1 spin up to a few percent, 
which freezes at low temperatures. The robust feature of 
the quadratic temperature dependence of the magnetic 
specific heat, and its scaling with the Weiss temperature 
indicate the existence of the Nambu-Goldstone mode of 
a gapless linearly dispersive type. Absence of either con- 
ventional magnetic order or bulk spin freezing suggests 
a novel symmetry breaking of the ground state. 

The polycrystalline samples of Nii_a;Zn2;Ga2S4 are 
synthesized by annealing of Ni, Zn, Ga and S elements in 
evacuated silica tube at 850 ^ 900 °C. In order to obtain 
a homogeneous mixture of Ni and Zn, we first ground 
the same molar amount of Ni and Zn powder, and re- 
peated the dilutions by adding equal molar amount of 
Ni powder each time until we obtain an appropriate Zn 
concentration. Powder x-ray diffraction analysis at room 
temperature on our polycrystalline samples indicates sin- 
gle phases with the same trigonal structure as NiGa2S4 
up to a; = 0.3. Our preliminary neutron diffraction mea- 
surements for X = Q, 0.1, 0.25 confirmed both the trigo- 
nal structure with P3ml symmetry and the substitution 
of Zn at the Ni site at the same concentration as the 
nominal one within an error of 1 Vo}'^ For x > 0.4, two 
phases coexistence of the trigonal and tetragonal phases 
is found because of a phase transition to the tetragonal 
phase of ZnGa2S4. In this letter, we focus on the physical 
properties of the trigonal phase region (0 < a: < 0.3). 

DC magnetization AI was measured between 1.8 and 
350 K under magnetic fields up to 7 T using a SQUID 
magnetometer. The specific heat Cp was measured by 
thermal relaxation method down to 0.35 K. In order to 
estimate the lattice contribution, Cl, we measured Cp of 
the isostructural nonmagnetic analogue ZnIn2S4 and get 
the thermal variation of Debye temperature Od (T) us- 
ing Debye equation. 6'd (T) of Nii_a;Zn2;Ga2S4 is then 
estimated by multiplying a scaling factor according to 



9d {T)(xM-''^Vc 



-1/2t/-1/3 



where Mq and Vq are molar mass 



and volume, respectively. Cl is thus estimated by con- 
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Fig. 1. Temperature dependence of the susceptibiUty X (T) for 
both FC (soUd) and ZFC (open) sequences under (a) B = 0.01 
T and (b) 7 T. The insets show (a) the x dependence of the 
freezing temperature Tf and Ax (T = 1.8 K) under B = 0.01 T, 
and (b) the effective moment p^g per Ni ion, the absolute value 
of the Weiss temperature |6w| and the frustration parameter /. 
The blue horizontal line gives p^ff = 2.84 expected for 5 = 1. 



verting back this scaled 9u (T) into specific heat. 

The temperature dependence of the susceptibihty 
X{T) = M (T) /B under fields B = 0.01 and 7 T is 
presented in Fig. 1. Under 0.01 T, samples of all concen- 
trations show the bifurcation between field cooled (FC) 
and zero-field cooled (ZFC) data below freezing temper- 
ature Tf (Fig. 1(a)). For the pure NiGa2S4, the freezing 
effect is so tiny that it is not the bulk spins that freeze, 
but nearly 300 ppm defect spins due to the imperfection 
of the samples including surface.^ The Zn substitution 
increases the difference between FC and ZFC data at 1.8 
K, Ax {T = 1.8 K), linearly with the Zn concentration x 
(inset of Fig. 1(a)). Under 7 T, no hysteresis between FC 
and ZFC data was observed down to 1.8 K (Fig. 1(b)). 

Above 150 K, the susceptibility data for all concentra- 
tion follow the Curie- Weiss law: xiT) = C/ {T-Bw)- 
The effective moment pos per Ni ion estimated from the 
Curie constant C is close to the value expected for 5 = 1 
(2.83 /is). This confirms that Zn substitutes Ni by the 
same amount as the nominal concentration. The Weiss 
temperature 9yf is found negative and AF. As in the in- 
sets of Fig. 1(a) and 1(b), the absolute values of both Tf 
and show a drastic change at a; ~ 0, indicating the 
strong sensitivity of the bulk disordered state of NiGa2S4 
to the impurity. Further substitution systematically de- 
creases both Tf and \9y/\ in a similar fashion. 

In order to estimate the amount of Ni spins affected 
by the Zn substitution, we performed a two-component 
analysis on the magnetization results. Here, the two com- 
ponents are the bulk spins most likely forming the spin 
disordered state similar to NiGa2S4, and "defect spins" 
induced by the Zn substitution. Comparing the results 
under 0.01 T (Fig. 1(a)) and 7 T (Fig. 1(b)), one notices 
that the low-T susceptibility becomes significantly field 
dependent with x, while NiGa2S4 exhibits nearly field- 
independent susceptibility. This suggests that the defect 
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Fig. 2. (a) Difference between the susceptibility under B = 0.01 
T and 7 T fitted by xo -I- yCs=\l (T - Q). The inset shows its 
fitting result, y vs. x. (b) Magnetization M vs. the field B fitted 
by eq. (1) in the text. The inset shows its fitting result, z vs. x. 



spins are much more weakly coupled than the bulk spins 
with |6'w| = 80 K, and can be fully polarized under 7 T. 
In order to confirm this, the difference between B = 0.01 
T and 7 T are fitted by Xo + yCs=i/ {T - 6) in the range 
of Tf < T < 20 K, where Cs=i is the Curie constant 
of 1.0 (emu K/mol) for S = 1 and a fitting parameter 
y gives the defect spin fraction. All curves are fitted as 
shown in Fig. 2(a). The Weiss temperature for defect 
spins, 9, was only about one half of Tf and much smaller 
than |^w|, confirming that the defect spins are weakly 
coupled. Xo is negligibly small ^ 10^^ emu/mol. The 
fraction of the defect spins, y, increases linearly with x. 

We also performed the two-component analysis on the 
field dependence of the magnetization M (B) . For the 
pure NiGa2S4, the field dependence data (0 < B < 7 T) 
exhibit a linear increase at 1.8 K. However, non-linear 
Brillouin function type component appears and increases 
with the Zn substitution (Fig. 2(b)) because of the weakly 
coupled defect spins. In order to estimate the amount 
of the defect spins, we performed the fitting using the 
following equation in the range of < S < 7 T, 

/ gfiBSB\ 



, dM 
MiB,x)^ — 



B + zgii^SBs 



B=7T 



(1) 



Here, we assume that the susceptibility for bulk spins is 
insensitive to field as in the pure NiGa2S4 and is given 
by at 7 T where the defect spins are polarized. Bg 
in the second term represents the Brillouin function for 
5* = 1 defect spins. The eq. (1) fits the experimental 
results well as shown in Fig 2 (b). The fraction of the 
defect spins, z, increases linearly with x. The above two 
types of the analyses consistently give the same order 
of the defect spin concentration that are proportional to 
X. This indicates that the nonmagnetic impurity sub- 
stitution induces S = I weakly coupled defect spins 
that most likely comes from Ni spins adjacent to Zn 
ions. Only a few percent of Ni site appears to induce 
the defect spins. Given the low temperature hysteresis 
Ax {T = 1.8 K) cx X, these results indicate that it is not 
the bulk spins but the defect spins that freeze. 
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Fig. 3. T dependence of (a) the total specific heat divided by 
temperature, Cp/T, (b) the magnetic specific heat devided by 
temperature, Cm/T. The inset shows the magnetic entropy 5 as 
a function of T/|0-vvl- The horizontal line indicates S = _Rln3. 




Fig. 4. T/|6»w| dependence of (a) Culdwl/T, (b) (Cm - 
Co)|6wP/T^, and (c) the derivative of the susceptibility dx/dT 
under B = 7 T. Co is defined for hD/Loks < T < l^wl- The 
vertical dashed lines indicate 4 % of |9wl a-nd Tq ~ 13 % of |9wl- 
The inset shows the x dependence of the coherent scale Lq and 
the stiffness constant D. Co = 0.0(2) mj/mole K for NiGa2S4 
indicates infinitely long Lq and its lower bound is plotted. 



The absence of bulk spin freezing as well as spin or- 
der was confirmed by specific heat measurements under 
B = T (Fig. 3(a)). The T dependence of the specific 
heat divided by T, Cp/T (T), shows only broad anoma- 
lies, indicating that conventional AF order is still sup- 
pressed in Nii_2;Zna;Ga2S4 and that the interacting spin 
remain disordered in the low-T limit. In order to probe 
the magnetic density of excited states, the T dependence 
of the magnetic part of the specific heat Cm {T) was esti- 
mated by subtracting the lattice part Cl{T) from Cp (T) 
(Fig. 3(b)). All X ^ compounds have a double-peak 
structure, similar to NiGa2S4 {x = 0). One is the broader 
peak centered around T ^ \0w\, and the other is the 
prominent rounded peak at the temperature close to Tf . 

Significantly, Cm/T for all x shows a T linear depen- 
dence at T — > 0. This is in sharp contrast with T indepen- 
dent Cm/T due to the local nature of spin fluctuations 
in canonical spin glasses. This indicates absence of spin 
glass state in the bulk, but the coherent propagation of a 
gapless and linearly dispersive mode in two dimensions. 

With a 2D gapless linearly dispersive mode with a co- 
herent propagation limited up to a length scale Lq at 
T = 0, the specific heat deviates from the low tempera- 
ture asymptotic form as 

Cm{T) _ Co , 3V3C(3) fakeTY 

R " i? + 27r \ hD J ' ^ ' 

for hD/LokB < T < {6^1, where Cq = 
- {V3tt/2) {a/Lof R and C (3) = f. 202.9^^2 jjgj.g ^ jg 
the lattice constant, and D the spin stiffness constant. 

For NiGa2S4, Cm shows quadratic T dependence with 
Co = 0.0(2) mJ/molc K and this indicates infinitely long 
Lq with the lower bound ~ 130 nm, which is far be- 
yond the two-spin correlation length ^ ~ 2.5 nm de- 
termined by the neutron scattering.^ With the substi- 
tution of Zn, however, Lq sharply decreases at a: '--^ 



(inset of Fig. 4). This shows the coherence of the pure 
NiGa2S4 is fragile and can be easily suppressed by a few 
percent impurity. On the other hand, the spin stiffness 
constant shows the continuous decrease with x (inset of 
Fig. 4). For ordinary AFMs that order at T |6'w|, the 
stiffness constant Dq can be estimated by the relation: 
Dl « (3%/3C(3)/47r) (afcB6'w/ft)Vln(25'-H l).^'i2 The 
observed D is between 0.65 and 0.85 km/sec, nearly three 
times smaller than the expected Do = 2.3 - 3.0 km/sec. 
This softening comes from the magnetic frustration, and 
indicates the spectral weight downshift. Further evidence 
for the downshift is found in the magnetic entropy 5* (in- 
set of Fig. 3(b)), obtained by integration of Cm/T. The 
observed intermediate temperature plateau of 5* prior to 
high-T saturation at i?ln3 amounts to 1/4 ~ 1/3 of 
i?ln3, indicating the highly degenerate low-T states. 

Interestingly, the observed systematic change of the 
specific heat and susceptibility as a function of x roughly 
scales with the Weiss temperature. When the system is 
dominated by the single scale |0w|i the molar entropy 
should take the form S = /(T/j^wD- Indeed, S pre- 
sented in the inset of Fig. 3(b) shows such a scaling as 
a function of T/\6w\, except one for NiGa2S4. Then, the 
specific heat should be given by T/\9w\ times a scal- 
ing function of T/|(?w|- In order to check this, we plot 
Cm\0w\/T vs. T/\ew\ in Fig. 4(a). Surprisingly, all the 
curves except the one for the pure NiGa2S4 overlaps on 
top of each other. Moreover, including NiGa2S4, the peak 
temperature of Cm\(^w\/T and the initial slope as a func- 
tion of T/\9w\ are constant, indicating that the peak 
temperature is proportional to |0w|i and T^ constant of 
Cm is proportional to |0w|~^- In order to further check 
this, (Cm-Co)|6'wP/T2 is plotted against ln(r/|6'w|) in 
Fig. 4(b). All the results including the one for NiGa2S4 
(with Co = 0.0 (2) mJ/mole K) starts to have the same 
constant value below the temperature of T/|6'w| ^ 0.04. 
While the lower T limit to define Co, hD/LokB corre- 
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spends to r/l^wl = 0.015(4) for the Zn substituted sam- 
ples, the eq. (2) roughly describes Cm down to the lowest 
T given by T/|6'w| ~ 0.004, much lower than the limit. 

On the other hand, the scaling behavior in the suscep- 
tibility is not so clear as in the specific heat results prob- 
ably because the susceptibility is more sensitive to defect 
spins at low temperatures. However, we note the three 
important characteristic temperatures that scale with 
|6'w|- (1) Below 4 % of |^w|) the susceptibility becomes 
constant, as indicated by the vanishing T derivative, 
dx/dT (Fig 4(c)). (2) The minimum of dx/dT appears at 
nearly the same scaled temperature of T/|6'w| ~ 0.2 (Fig 
4(c)). (3) Tf scales with \9w\, so the frustration param- 
eter / = \6w\/T{ keeps almost a constant value around 
10, indicating the significant magnetic frustration even 
with the Zn substitution (inset of Fig. 1(b)). 

Despite the above scaling behavior with the Zn substi- 
tution, only the pure NiGa2S4 exhibits the following two 
noticeable anomalies at the crossover temperature to the 
low-T asymptotic behavior. (1) Not following a scaling 
function for the Zn substituted systems, only NiGa2S4 
exhibits a broad kink in (Cm - Co)\Ow\^/T^ at Tq ~ 13 
% of |6'w| — 10 K, below which nearly dependent Cm 
appears. (2) Only for NiGa2S4, dx/dT crosses zero at 
the same characteristic temperature Tq , below which the 
nearly constant x appears. Interestingly, this low tem- 
perature crossover behavior in the pure NiGa2S4 can be 
easily suppressed only by 1 % Zn, similar to the rapid 
loss of the coherence scale Lq. 

The Zn substitution substantially but never com- 
pletely suppresses Lq and dependent specific heat. 
The observed robust scaling of the low temperature 
behavior, especially the form of the specific heat, 
strongly suggests the existence of the Nambii-Goldstone 
mode, which has a gapless and linearly dispersive charac- 
ter and a coefficient scaling with j^wl"^-"^^ Generally, 
a Nambu-Goldstone mode appears with a broken sym- 
metry in comparison with its high-T phase. Therefore, 
in our case, the absence of the evidence for conventional 
long-range order in thermodynamics and neutron results 
points to a novel magnetic order in two dimensions. 

One candidate for such a unconventional long-range 
order without long-range two-spin correlation is a spin 
nematic phase. ^'^ This phase may be classified as a spin 
liquid because the static site average of spin is zero, while 
its magnetic quadrupole correlation is long-ranged. If the 
transition occurs at T = Tq, spin fluctuations should 
show a critical slowing down, which may well lead to the 
freezing of impurity spins at Tf close to Tq. Recently, the 
mean-field type calculations by Tsunetsugu and Arikawa 
have shown that a spin nematic order can be indeed sta- 
ble on a 2D triangular lattice, and generates a Nambu- 
Goldstone mode of the gapless linearly dispersive type.^^ 
In our case, the strong two-dimensionality of the spin in- 
teractions may suppress the three-dimensional character 
of the order and instead promote a rapid development of 
the correlation at a crossover temperature Tq, resulting 
in the broad feature of the specific heat near Tf. Fur- 
thermore, the orthogonal spin nematic order with three 
sublattices has no geometrical frustration on the triangu- 



lar lattice,^'* and may be robust against the substitution 
of a sizable amount of nonmagnetic impurities. 

Another candidate is a Kosterlitz-Thouless (KT) type 
phase driven by two-valued vortices with the transition 
around Tf.^^ Slow dynamics associated with bound vor- 
tex like defects may cause non-ergodicity below a KT 
type transition into a critical state with a finite ^. 

Other possibilities include the critical behavior toward 
an AF transition below T = 0.35 K.^^ However, it is not 
yet clear whether such a low-T AF order with geometrical 
frustration could be robust against nonmagnetic impuri- 
ties of the concentration up to 30 %. 

Finally, a quantum spin liquid phase is also an interest- 
ing candidate. For triangular AFMs, Imada et al. discuss 
a gapless spin liquid phase through the theoretical cal- 
culations based on the Hubbard model. However, it is 
not yet clear whether the formation of such a liquid state 
may lead to the T^ dependence of the specific heat. 

To conclude, the Zn substitution for Ni in NiGa2S4 
leads to the abrupt loss of the coherence by the substitu- 
tion less than 1 %. However, the loss is not complete and 
the robust low temperature behavior, i.e. the constant 
susceptibility, the T^ dependent specific heat, and their 
scaling behavior with the Weiss temperature indicate the 
existence of Nambu-Goldstone mode of the gapless lin- 
early dispersive type. The absence of either conventional 
magnetic order or bulk spin freezing suggests that the 
ground state has a novel symmetry breaking. 
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